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ABSTRACT 


An analytical and test program was conducted in order to establish the feasi- 
bility of a multichamber pulsed-capillary viscometer. The initial design 
incorporated a piston, ram, and seals which produced measured pulses up to 
30,000 psi (206,800 NT/M 3 ) in the closed chamber system. Pressure pulses 
from one to ten milliseconds were investigated in a system volume of 1 in 3 
(1.63 x 10~ 5 M 3 ). 

Four test fluids: a MIL-L-7808, a 5P4E polyphenyl ether, a MIL-L- 23699A, and 
a synthetic hydrocarbon were examined in the test pressure assembly. The 
pressure-viscosity coefficient and viscosity delay time were determined for 
the MIL-L-7808 lubricant tested. 
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INTRODUCTION 


The operation of rolling element bearings and other high speed concentrated 
contact elements such as cams, gears, and rolling contact drives involves sub- 
jecting lubricants to extremely high pressures of the order of hundreds of 
thousands of psi and high shear rates for very short duration times. Many 
investigations have been performed to attempt to understand the behavior of 
lubricants under these conditions* Such investigations include: 

1. High pressure viscometry which provides measurements of lubricant 
viscosity under high pressures and nearly isothermal conditions and 
at low shear rates and long times* 

2. Rolling- sliding contact measurements which involve obtaining 
traction and film thickness data with a highly loaded rolling disc 
machine or other type of rolling contact apparatus* Here the lub- 
ricant is subjected to high pressures as well as high shear rates, 
short times and often high temperature rises* 

This second class of measurements does simulate the conditions undergone in 
concentrated contact elements* However, they simulate all of the conditions 
simultaneously (high pressures, high temperatures, high shear rates and short 
times) with a geometry that is somewhat uncertain in that the film thickness 
is determined by the lubrication process itself and must be measured* 

The results of these rolling contact measurements indicate at the high speeds 
encountered in a concentrated contact element, various phenomena occur that are 
not explainable by high shear rate and thermal effects alone* In particular, 
variations of measured traction and film thickness with rolling speeds can not 
be accounted for by long time equilibrium Newtonian theory nor can they be de- 
scribed by simple shear rate dependent lubricant models* It has been hypothe- 
sized that short time effects could be responsible for some of these observa- 
tions* In particular, a lubricant moving at a thousand inches per second will 
pass through a 10 mil contract zone wit hr a ‘ resldenceu time ' 6f 1.0- microseconds . 
During this passage through the contact zone a particle of fluid will experi- 
ence a change in pressure from ambient to hundreds of thousands of psi and 
back to ambient* It has been hypothesized that the viscosity which increases 
radically with pressure during long times does not increase quite as radically' 
during short times* In other words, there is a delay time associated with the 
response of the viscosity to a rapid change in pressure* This viscosity delay 
phenomenon is the subject of the present investigation* 

In order to study viscosity delay phenomenon, a simple apparatus consisting of 
a pulsed capillary viscometer has been developed* A pressure pulse is trans- 
mitted to a reservoir full of fluid and a small amount of fluid is forced 
through a capillary to another closed reservoir full of fluid* By measuring 
the pressure time characteristics of both reservoirs it is possible to infer 
the viscous characteristics of the fluid* Pulse times can readily be achieved 
to the order of a few milliseconds and residence times in the capillary tube 
can be as small as a few microseconds* The capillary tube can be built to 



have a relatively short length and small but controlled radius hence pro- 
viding a controlled geometry. It was hoped that this device could provide a 
simplified means of measuring viscosities under moderately high pressures and 
short times but at lower shear rates with smaller thermal effects than en- 
countered in rolling contact apparatus. Hence, it is that the pulse type 
viscometer will have the advantage over the disc machine type of apparatus 
in that it provides a controlled geometry and at least partially isolates the 
short time effects from high shear rate effects in a relatively simple apparatu 

The following sections describe the concept of the pulse viscometer, the 
analyses required to relate pressure pulse measurements to viscosity, viscosity 
pressure and viscosity- time characteristics, a description of experiments that 
have been performed together with results of a viscometric analysis of the 
test data. 
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THE PULSED- CAPILLARY VISCOMETER 


Conceptual Discussion 

High pressure viscometers (1) which incorporate capillary tubing as flow re- 
strictors have been in use since 1966. Viscometers incorporating capillary 
tubes in the past have been of the slow flow variety. Fluid under high 
static pressure is induced to flow between two chambers via capillary tubing 
under a relatively small differential pressure. 

The present design incorporates two or more chambers which are initially un- 
pressureized and connected by capillary tubing to the largest volume chamber. 
The design concept is shown schematically in Figure 1. The design is unique 
in that flow through the capillary tube is produced when the primary chamber 
is pressurized momentarily by impacting the enclosure piston. Each chamber 
of the viscometer contains an active pressure transducer element which provides 
a readout of the pressure-time history within each cavity of the high-pressure 
system. An oscilloscope display of a typical pair of pressure pulses from a 
single impact is shown schematically in Figure 2. The impact is produced in 
such a way as to provide a half-sine wave pressure-time pulse in the primary 
chamber of the viscometer. The pressure buildup as a function of time in the 
secondary chamber is characteristic of the rheological properties of the con- 
tained fluid in the closed system. Correct sizing of volume V 2 and the capil- 
lary tubing leading to the primary chamber provides the proper readout sensi- 
tivity for the viscosity of any fluid chosen for investigation. 

Capillary tubes which are too large will allow the* secondary pressure (P 2 ) to 
follow the primary pressure (p-^) and no estimate of fluid test viscosities 
could be made. Conversely, capillaries which are too small will not allow 
sufficient fluid to enter the secondary chamber during an applied pressure 
pulse and no signal would be obtained from the secondary chamber. Proper siz- 
ing of the capillary tubes is required for the desired output sensitivity and 
can be estimated from the approximate test fluid viscosity under investigation. 

The relative size of the test volumes, and V2> as shown in Figure 1 is also 
important when considering the desired range of pressure ratios P 2 /PI (see 
Figure 2). If V 2 is too large, the filling time will be extremely long and 
will not provide sufficient pressure buildup in the secondary chamber. A small 
volume V 2 would create a condition similar to an extremely large capillary in 
that the pressure ratio p^/p^ would be one and the instrument would provide no 
usable information. 

A consideration of the base viscosity of the test fluids is necessary to deter- 
mine the proper secondary chamber and capillary sizes. Highly viscous fluids 
will reduce pressure buildups in the secondary chamber as a result of reduced 
flow. On the other hand, fluids with low viscosities might not be flow resis- 
tant enough when passing through improperly sized capillaries and would provide 
little or no viscosity readout by allowing (p ) to follow (p^) too closely. 

Properly sized chambers and capillaries for the expected fluids to be tested 
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will yield a peak pressure in the secondary chamber which is one-half or less 
of the peak value attained from a pulse applied to the primary chamber . The 
analysis described in detail in the following section prescribes the parametric 
interdependence for the high-pressure pulsed-capillary viscometer and gives an 
insight into how one can properly size chamber components for correct sensi- 
tivity. 

Test Assembly and Data 

Although the test approach used in measuring viscosities with pulsing techniques 
is new, the test apparatus is a combined assembly of commercially available and 
new machined components. The nature of the physical requirements dictated spec- 
ific design features which were chosed for testing and are reviewed in this 
section. 


Primary Pressure Chamber 

The primary chamber (Volume V^ in Figure 1) is contained in a commercially 
available cylinder capable of withstanding static pressures up to 200,000 
psi 0 The test cylinder was: constructed of two prestressed walls and has 
a six inch outside diameter and a five-eighths inch bore. The piston ram 
during operation drives a floating seal and plunger into the bore of the 
chamber. The standard seal used in the plunger was replaced with an elas- 
tomeric Oaring in order to ensure a positive seal with a minimum of sliding 
friction. Visual evidence of sealing (at least to the pulsed 30,000 psi 
pressure level) was shown by the return of the piston to its original 
height after every impact. Perfect sealing of the primary chamber was 
not necessary for operation since the primary pressure buildup was moni- 
tored during pulses. 


Secondary Chambers 


The theory of the pulsed viscometer concept (to be described later) provides 
a dimensionless parameter g which prescribes appropriate dimensions of the 
capillary and secondary chamber combinations. An appropriate combination 
of parameters requires that (3 < 1 for proper sensitivity where 


P = tt R 4 kT 
8 |i LV 2 


( 1 ) 


The diameter of the capillary tubes used in the high-pressure tests was 
determined experimentally from a set of low pressure flow tests. The 
integrally mounted capillary and associated secondary chamber was cali- 
brated by forcing water through the tubing. Measurement of the flow rates 
of distilled water through the known length of capillary tubing provided 
an accurate determination of the radius of the capillary tubing. In 
addition, volumetric determination of the test secondary chambers from 
machined dimensions provided the remaining information required to evaluate 
the geometric portion of g defined below: 


geometry factor = 


TT 

8 L V 2 
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for each test setup used in obtaining the oscilloscope pressure data 
presented in this report. The calibration for each test capillary indi- 
cated 

Test Capillary Designation 
A 
B 
C 
D 
E 

Test capillary tubes were produced from sectioning standard stainless steel 
hypodermic tubing. Lengths varied from .100 to .400 inches and the inside 
diameter from 4 to 8 mils. The securing of the tubing in line between 
chambers was accomplished with heat cured epoxy and worked well over the 
pressure range tested in the present program. 

The primary pressure chamber contained approximately one cubic inch of test 
lubricant. A determination of the precise volume was not important as 
long as the primary chamber volume was at least twenty times greater than 
the secondary volume. The primary chamber in the present study was 25 to 
100 times that of the secondary volume depending upon the initial position 
of the piston ram. 

Pressure Gages 

Two types of pressure gages were tested during this program. Both were 
made of manganin sensing elements, one in the form of a wire wound open 
coil and the other in the form of a thin 200 micro-inch foil. Although 
each type of sensor was available commercially, each type had to be modi- 
fied slightly for mounting on the secondary chamber follower and seal. 

Figure 3 shows a mounted coil of the wire type sensor ready for operation. 

A second coil identical to the one shown was mounted externally to the 
pressure chamber to form a half-bridge in a standard Wheatstone arrange- 
ment. When mounted in the Wheatstone bridge, each sensor provided a 
100 millivolt output at 10,000 psi pressure. The resolution of the oscil- 
loscope pressure- time trace, as a result, was limited to + 25 psi. Pulse 
periods from two to ten milliseconds were determined with an accuracy of 
5 percent or better from the oscilloscope tracings. 

Foil gages made of manganin were also evaluated for pulsed pressure deter- 
mination. The foil was an open etched pattern 200 microinches thick and 
had an active grid area of approximately 1/8 inch square when bonded to 
the chamber follower as seen sectionally in Figure 4. Again the output 
sensitivity of the element was adjusted to be the same as the wire type 
gage and was statically calibrated with the aid of a hand pump at 10,000 
psi in a closed pressure chamber. A standard Bourdon tube pressure gage 
with 10,000 psi scale aad a + 1 percent accuracy was used for calibrating 
the pressure detection bridge. 


Non-Dimensional Geometry Factor 
2.027 x 10' 8 
3.544 x 10' 8 
1.418 x lO -9 
1.229 x 10~ 9 
1.893 x 10“ 9 
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Of the two types of active pressure elements used, the foil gage was pre- 
ferred for its small size and minimum volumetric displacement* It is 
possible that a modified version of the foil pressure sensor could be 
incorporated within a capillary channel if analytical evidence was pro- 
duced to require such a measurement. 

Instrumentation 


A schematic of the instrumentation setup is presented in Figure 5* The 
simplicity of the pulsed pressure viscometer is apparent from the diagram* 
The essential components are pressure transducer elements^ Wheatstone 
bridges s and a multi-beam oscilloscope with camera attached* Photographs 
of Figures 6 and 7 show the bridge instrumentation and the high-pressure 
cylinder. Each bridge contains its own zero balance control and amplifier 
with span adjust for calibrated voltage output. The input connections of 
the bridges used allowed hookup for either two or four external bridge arm 
connections. The signal to noise ratio for the smallest pressures measured 
was 20:1 whereas the ratio went to a 1000:1 for the higher pressure signals 
detected* The frequency response of the bridge system was limited to the 
zero to twenty thousand cycles per second range in the present facility. 

The necessary bridge frequency response for the shortest pulses induced in 
the present study was one thousand cycles per second. A full list of 
instruments used and their model numbers are included in Appendix 1 for 
further reference . 


Test Data 


The output data from the present test program is contained in over one 
hundred and fifty dual trace pressure-ver sus- time oscilloscope photographs* 
Pressure- time traces for each of four fluids were recorded. The evaluated 


fluids were 

Fluid 

MIL-L-7808 
Polyphenyl Ether 
MIL-L-23699A 
Synthetic Hydrocarbon 

and some of their physical proper tie 
performed prior to the present study 


Description 
Humble #2389 
Monsanto OS-124 
Mobil Jet II 
XRM 214B 

s as determined from investigations 
are included in Appendix 2. 


In addition to determining chamber pressures during test impact^ checks 
were made of the test fluid compressive heating with the aid of small 
thermocouples. Two five-mil diameter thermocouples were joined and flatten- 
ed to less than a 2-mil thickness. Stepped thermal changes applied to the 
produced thermocouples revealed a response time of approxima tely three 
milliseconds. Using the miniature temperature sensors^ the bulk compres- 
sive heating within the secondary and primary chambers were found to be 
less than 4°F for applied pressure pulses of approximately 20 9 000 psi 9 
lasting from 5 to 6 milliseconds. 



Viscosity data are inferred analytically from the pressure- time traces 

obtained from the primary and secondary chambers after impact „ Typical 

pairs of traces obtained for each lubricant tested are displayed in 

Figures 8 through 31* The pressure (p, ) of the impact load is found in 

each run from the upper trace, depicted in each photographo Since the 

scale factors vary from picture to picture, the peak pressure in the 

primary chamber (p^ ) and the peak pressure in the secondary chamber 

(p 0 ) are noted at ?he side of each run shown*, 

max 



PULSE VISCOMETER ANALYST S 


Pulse viscometry measurements consist of input temperatures, times, pressures, 
and geometric variables including reservoir volumes and capillary lengths and 
diameters. In order to convert these measurements into fluid properties, it 
is essential that analyses be performed which relate fluid properties (mainly 
the viscosity and the viscosity pressure, temperature, and time characteristics) 
to these pressure measurements. 

The analyses to be presented here are divided into two basic parts. The first 
of these represents a development of simplified relationships assuming highly 
idealized viscosity-pressure- tempera ture relationships and no time effects 
which has been used for purposes of designing the test apparatus and extract- 
ing basic preliminary information about the relationship and sensitivity be- 
tween measured pulses and long-time viscosity characteristics. The second set 
of analyses is somewhat more comprehensive in that it allows for a viscosity 
delay, it allows for an arbitrary input viscosity-pressure-temperature relation- 
ship (for equilibrium viscosity), and allows for other variable properties such 
as bulk modulus variations with pressure,, 

The first set of analyses thus provides qualitative trends and ranges of rele- 
vance for the various parameters as well as identifying relevant parameters. 

The second set of analyses provides a more sophisticated tool for extraction 
of viscosity and viscosity delay characteristics. 

Simplified Analyses for Predicting Long-Time Pulse Measurements 

A schematic showing the two reservoirs, noted as reservoirs 1 and 2 for inlet 
and exit ends of the capillary respectively, are shown in Figure 22. A pulse 
is imparted to reservoir 1 which acts as a spring, providing an approximately" 
sinusoidal pressure pulse which forces fluid through the capillary into reser- 
voir 2. In this section the flow through the capillary and the resulting 
pressure profile, P 2 (t), will be analyzed first for an isoviscous fluid and 
then extended to include viscosity-pressure and viscosity-temperature effects. 

Isoviscous Fluid 

The volumetric flow rate, Q, of an isoviscous fluid undergoing laminar 

flow through a capillary is given by 

0 _ TTR 4 (p i " P2 ) (2) 

^ 8\1 L 

If the density pressure relationship is of the form 

P - P Q (1 + £) (3) 
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where k is the bulk modulus of the fluid and p/k < 1 ( .2 or 

less) then the pressure, is related to the flow, Q, by 





v 


2 


If the input pressure (p^) 


is of the form 


p i 


af ( t / t 0 ) * t<t o 

0 , t > t Q 


( 4 ) 


(5) 


where t denotes the duration time of the impulseand f is a dimensionless 
input p§essur£ function then we may combine Equations (1), (2), and (5) to 
obtain the relationship 



4 

ttR k 
8yLV 2 


[af(t/t o > - 


P 2 ] 


0 < t < t 

O 


( 6 ) 


Now introduce dimensionless variables 


P 2 = p 2 /a and t = t/t o 

and the parameter 
4 

ttR kt 

g = ° 

8pLV 2 

to obtain the dimensionless equation 

0 < t < 1 
t > 1 

P 2 (0) = 0 

The above first order equation has the solution 


d?2 

dt 


8 { [f (t) - 


- 8p, 


p 2 ]} » 


C7) 



-8 (t-t 1 ) 
e 


f (x 0 ) dr't 9 




8t’ 


f ( t 5 )dx 


> 


T < 1 


’ T>;> 1 


( 8 ) 
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If mechanical energy, E, is imparted to a system of mass, M (M would be 
the mass of the piston if the projectile did not stick to the piston 
during the pulse time or the sum of the masses of the piston + projectile 
if sticking did occur), then the pulse pressure would be 



where A denotes the cross sectional area of the piston* The above formula 
neglects sonic effects, assumes no leakage, and is strictly valid for 
V » V 2 * Equation (9) was arrived at by considering a simple spring and 
mass system . An initial velocity v is imparted to the system given by 
v Q = /2E/M. If the piston displaces a distance x, the volume, displaced 
would be AV]_ - Ax and for small displacement ~AV]_/V^ » Ap/p 0 = p/k* The 
right hand side of the above equation is a consequence of Equation (3). The 
spring stiffness S would then be given by S = Ap/x = RA^/Vi- Equation (9) 
thus represents the simple spring and mass solution for a system having mass 
M, stiffness S given as initial velocity v given by the above expressions. 
The amplitude of the pressure pulse is 


\ f 2kE 

IK 


( 10 ) 


and the period is 


MV 


t = TT 
O 


kA 


(id 


hence 

f (t ) = sin TTT 

The solution to Equation ( 8 ) with f (t) = sinTTT is 


— 2 O sinrrT " ttcosttt + rre ^ T ) 0 < t < 1 


( 12 ) 


? 2 ■ ( 


TT + 


~~ 777 2 


(13) 


T > 1 


The value of t for which p 0 is maximum is given by the solution to 


e(e 


max 


COSTTT ) 

max 


TTsinTTr 

max 


(14) 



Values of T max calculated from Equation (14) along with corresponding 
values of p max are shown in Figure 23 0 It can be seen that for greatest 
sensitivity of measured pressure, p^> to viscosity, the parameter, (3, should 
be considerably less than one, " 

The Effect of Viscosity-Pressure Dependence 

If the viscosity varies with pressure in accordance with the relationship 
P = P Q g(p) 

then Equation (1) assumes the form 


Q 



(15) 


and Equations (4), (5), and (15) may be solved simultaneously* If we 
assume a viscosity-pressure relationship of the form 


g(p) *= 


(16) 


where a is the conventional pressure coefficient of viscosity then Equa- 
tion (15) becomes 


Q 




- e 



(17) 


and the differential equation for p^ is as follows 


d P. 

z 

dT 


-cvp 9 -ad: (t) 

(e - e ) , 


a 


- 


a 

~ (e -1) 

a 


0 < T < 1 


T > 1 


(18) 


P 2 (0) - 0 

where the parameter a is defined by a = aa . 

The above equation may be put into linear form with the transformation 
“P ? 

0 = e (19) 

as follows 
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-at (t) 

0(1 - 0e ) , 


^ pci - 0) 


0 < T < 1 


T > 1 


Equation (20) has been solved numerically for 0 for f(T) = sinTTT. Maxi- 
mum values of are plotted against |3 for various values of a in Figure 
24. 

An asymptotic expression has also been obtained for the solution to 
Equation (20) in the interval 0 < t < 1 based upon large values of the 
parameters a and ®tt/| 3. The resulting expression for 0 is 


0 . i + ex --- [i 


-QttTT -QTT(I-t) 

5 + (1 + 0)e ] 


and the maximum value of 0 occurs at 


An I 1 + 


Hence, for the major portion of the pulse duration 0 will be linear with t. 
If (3 is small and a is large then p^ will approach the form 


12p Q aLV 2 


0 < t < T 


hence, P 2 will depend upon the product of the base viscosity (|a Q ) and the 
pressure coefficient of viscosity (a) which is similar to the dependence 
predicted for film thicknesses in ela stohydrodynamic contacts. 

Inclusion of Thermal Effects 

The nonintegra ted form of Equation (15) that is valid when \i varies with 

axial posit ion, z, but does not vary with radial position,r, is obtained 

simply by replacing (p -p )/L appearing in Equation (2) by -dp/dz 

1 2 


Hence , 


4 

ttR d^ 
8|_i dz 


and the average velocity, u, is given by 


e _1L. 


RL ,&>.) 

4n Mz' 


12 


The volumetric rate of power dissipation due to viscous heating (neglect- 
ing radial dependence on p) is given by 



$ .ail. 2 r 4 ( ds_ 2 

$ v = * ( I7 } s 4^ <3P 


§, 


Hence, the heat generated per unit length L is given by 


r R .42 

/ r $ v dr s ^ (Jr) 


2 TT 


a 


8p dz' 


( 24 ) 


If the average rate of heat generated per unit length is equated to the 
rate of removal by convection then 

„ dT * 

P c 0 Tz ‘ § L 

and when u and are replaced by Equations (23) and (24) the resulting 
expression is 


pc 


dT 

dz 


dp 

dz 


which may be integrated to yield the equation 


T 



(25) 


where 1^ denotes the temperature in reservoir 1. 

Equation (25) thus indicates that subject to the limitations imposed by 
radial averaging and the assumed modes of heat generation and removal 
the temperature rise depends only on the heat capacity and the local 
pressure drop and is otherwise independent of viscosity and geometry. 

If the viscosity varies with pressure and temperature in accordance with 
the relationship 

H = |i g (p»T) (26) 

o 

then the effect of temperature variations may be accounted for by replac- 
ing Equation (15) with 


Q - 


ttR 

8u L 


dr 

* Pi-p 

8 (P, T l + — 


In the special case 


e 


ap-6(T-T 1 ) 


(27) 


(28) 
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The denominator in the integrand in Equation (27) assumes the form 


g (Pi 


P£P 

Pc 


fa 

pc 


<« + 7^ )p 


•) 




and the differential equation for pressure analogous to Equation 
becomes 


(18) 



dT 


— 6 ^ -(^ + *) p 2 + * f ( T ). e -« f ^)], o < T < 

(a + 1 >) 


( 29 ) 

T 

max 


where the parameter 1> is defined as 



(30) 


Hence, with the use of the viscosity temperature model given by Equation 
(28) it is possible to represent the maximum pressure ratio (P^ax^ SS 
functions of the three parameters P,c?, and Curves showing the vari- 
ation of P 2 max 0 at various values of a and 6 are given in Figures 

25 through 27, 


It can be seen from Figures 25 through 27 that in all cases p 2 max becomes 

insensitive to {3 as (3 becomes large* In that (3 is the sole parameter 
reflecting the viscosity at ambient pressure it is essential that sig- 
nificant sensitivity with respect to (3 prevail, hence, based upon these 
figures one should design the apparatus with values^of (3 somewhat less 
than one* The practical range of parameters o' and 6 for present test 
conditions with present lubricants should be between 0 and 3 for both 
parameters* Bearing this in mind, it is possible to use Figures 25 
through 27 to design the test apparatus so that necessary ranges of op- 
eration are encountered where (3 is sufficiently small to provide sensi- 
tivity to viscosity yet large enough to provide measurable values of p^. 

In order to more rigorously extract viscosities and viscosity delay in- 
formation, one should have an analysis that incorporates the best avail- 
able viscosity temperature and viscosity pressure data and bulk modulus 
data along with an inclusion of the effects of a viscosity delay* In the 
following section such an analysis is developed* 


Refined Analysis Including Viscosity Delay Effects 


Viscosity Delay Mechanism 

A first order lower viscosity delay mechanism has been used to quantify 
and extract time effects from the apparatus* 


Suppose at time t = 0 a fluid previously existing under pressure, p , is 

a 
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suddenly exposed to pressure, p^ 0 If the long-time equilibrium viscosity 
pressure relationship is denoted by p (p), then the instantaneous fluid 
viscosity is assumed to rise at a rate which is proportional to the de- 
parture from equilibrium in accordance with the first order rate equation 


, . P " P e ( P b > 

&L s - 

dt t, 

d 


(31) 


where the constant t^ has dimensions of time and will be referred to as 
the fluid delay time. 

The solution to Equation (31) subject to the initial condition \i = [X (p ) 
at t * 0 is 6 3 


-t/t 

M- = ^ e (P b ) " Cp e (P b ) " P e ^ P a^ e ( 32 ) 

of the variation of j j. with time in accordance with Equation (32) is shown 
schematically in Figure 28. 

Since viscosity is a function of temperature as well as pressure, it is 
necessary to hypothesize a model for the time response due to a sudden 
change in temperature. For the present time, it will be assumed that 
viscosity responds instantaneously to temperature and our attention will 
be confined primarily to the delay associated with a rapid change in 
pressure. 


We will confine our attention to an equilibrium viscosity-pressure-temp- 
erature relationship of the form 


M- e = G e (p) F (T) (33) 

where is the equilibrium viscosity at ambient pressure and reference 
temperature T Q and G e represents the pressure portion of the equilibrium 
(long-time) viscosity function and F(T) reflects the temperature varia- 
tion in viscosity. 

The instantaneous viscosity is given by 

H = G(p) F(T) (34) 

where G (p) is assumed to follow a first order delay law analogous to 
Equation (31) 


dG 

dt 


(G - G e ) 


and F(T) is assumed independent of time. 


(35) 
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Equations for Predicting Effects of Viscosity Delay on Pressure Measurements 

Equation (35) provides the necessary relationship for predicting viscosity 
delay effects in reservoir 1 where the fluid experiences a pressure pulse, 

P-i (t), applied over a duration of time, t Q . 

Hence, in reservoir 1 and at the entrance to the capillary (z = 0) the 
viscosity delay equation is 



(G 1 ~ G e (p l } , 


x = 0 
0 < t < t 

o 


where the subscript 1 denotes reservoir 1. 


(36) 


Equation (36) applies to jj. as well as G here since compressive heating is 
small and there is no significant shear heating in the reservoirs. Hence, 
the reservoirs should be very nearly isothermal at temperature 

In the capillary, however, a particle of fluid passes through a pressure 
field and the time derivative, d/dt, is replaced by the "total derivative" 
u d/dx. If one again replaces the radially-varying velocity, u, with the 
average value, u, as defined in Equation (23), the delay equation in the 
capillary becomes 


u 


6G 

6z 


G - G (p) 


0 < z < L 


(37) 


The flow rate Q is again related to the pressure buildup in reservoir 2 
$2 by Equation ( 4 ) and is applicable to a variable as well as constant 
bulk modulus. When the flow rate given by Equation ( 4 ) is equated to 
that given by Equation (22) and the resulting equation is solved for the 
pressure gradient, the following expression is obtained 


6z 


8^V 2 p 2 

________ 

ttR k(p^) 


0 < z < L, t > 0 (38) 


The average velocity 
of p in, accordance 


, u , 

wi th 


appearing in Equation (37) is expressed in terms 
Equation ( 4 ) to obtain the differential equation 


6r ttr 2 k (P?) 

^2 -6z S “ V^~ t d (p) ‘ G e > 0 < z < L (39) 

The temperature profile in the capillary resulting from viscous heating 
is still related to pressure by Equation (25), hence, the viscosity p 
appearing in Equation (38) may be expressed in terms of G and p as 

H s Hj. G(p) - pc (T T^)] (40) 


16 



The reference temperature, T , has now been taken to be the inlet reser- 
voir temperature, T^. 

At each instant of time. Equations (38) through (40) constitute a non- 
linear second order system with two dependent variables (G and p) and 
one independent variable z 0 The equations also contain one unknown para- 
meter p^. 

The pressure, p^, is a measured input variable which is either interpol- 
ated from a table or approximated with a sinusoidal input function* 

Other input quantities are geometric parameters R, L, and V^ 9 viscosity 
functions G e (p) and F(X), bulk modulus function k(p) and delay time func- 
tion t^(p). The principal output quantity will be the pressure p^- 

Equation (36) may be solved numerically by the fourth order Runge-Kutta 
process subject to the boundary condition at t = 0, G - 1 to obtain G^ 
as a function of t a 

There are three constraints on the second order system given by Equations 
(38) through (40). The first two of these are the conditions at x “ 0 

G = G^ 9 p - p^ at x « 0 (41) 

The third condition at x = L is 

p ~ P 2 a t x = L (42) 

which determines 

Equations (38) through (40) are solved numerically at each time step with 
the fourth order Runge-Kutta process. The value of p^ is varied system- 
atically in a nesting procedure until Equation (42) is satisfied. The 
value of p^ at the next time increment is obtained with the central dif- 
ference formula 


P 2 (t n + 2At P 2 (t n } " P 2 (t n-l ) 

A time sharing computer program is available for solving the above equa- 
tions to any reasonable degree of accuracy for any given set of input 
functions and parameters. 


Results of Viscosity 


Computa tions 


In order to qualitatively illustrate the effect of viscosity delay and 
to estimate the sensitivity of the measured pressure, p^, to a viscosity 
delay, solutions to the system of equations were carried out for the 
isothermal case, for a fluid having a constant pressure coefficient vis- 
cosity constant bulk modulus, and constant delay time (t^). The viscos- 
ity functions thus become 


(T) - 

(p) = e' 
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Also, a sinusoidal input was used as in the previous parametric studies. 

Under these conditions non-dimensionalization shows (t) to be a func- 
tion of four parameters. These are the parameters 3 and a used previous- 
ly and two additional parameters and X defined as 


- t/t 
d o 


(A3) 


and 


X 


t a 
TT d 

8 n- 



( 44 ) 


The above two parameters are the effective delay parameter in the reser- 
voir 1 (t,) and the delay parameter in the capillary (X), The ratio of 
d 

the two parameters 



8 (L) 

TT at 

o 


2 


for a typical case corresponding to a pressure amplitude of a “ 20 ksi, a 
capillary length of L = 0,3 inch, a radius of R = 2 mils, a duration 
time of t = 3 milliseconds and a viscosity of \i = 25 cp is calculated 

to be 3,5 x 10 3 . The fact that this ratio is extremely small compared 
with the duration time of the pressure pulse; hence, for the parameter 
ranges of current interest, one of two extreme cases can be expected to 
occur. These are: 


1. The delay in the capillary is of the order one or less, hence, 
there is no appreciable delay effect in the reservoir (t^-* 0) „ 

2o There is appreciable delay in the reservoir, hence, complete 
delay occurs in the capillary (the fluid viscosity in the cap- 
illary thereby being equal to the inlet viscosity) which would 
correspond to X -» oo 8 

Cases 1 and 2 are shown for typical values of p and a in Figures 29 and 
30 respectively. It can be seen that in either case there will be con- 
siderable sensitivity to the delay parameters. 

In general, delay times are expected to be small for the fluids, pressures 
and temperatures presently under consideration in comparison with the 
pulse time, t G . Case 1 is thus anticipated to prevail with no signifi- 
cant delay effects in reservoir 1 and with delay effects in the capillary 
governed by the magnitude of the parameter X, 


The dwell time of a particle of fluid in the capillary can be regulated by 
regulating the length of the capillary as evidenced by the (R/L) depen- 
dence of X. 



VISCOMETRIC ANALYSIS OF TEST DATA FOR A.iMIL-L- 7808 . FLUID 


A summary of pertinent data for twelve different sets of pressure measurements 
for a MIL-L-7808 fluid is shown in the first four columns of Table 1 . The 
first column of Table 1 contains the run number; the second column contains 
the amplitude of the inlet pressure^ a 2 in thousands of psi; the third column 
contains the measured duration time* t QS , in milliseconds; and the fourth column 
contains the amplitude of the pressure in the discharge reservoir P2 m ax S) 
thousands of psi* 

The pressure pulses* p^(t)* were very nearly sinusoidal; hence & the amplitudes 
were sufficient for describing the entire curve along with the period* t 0 « The 
remaining columns of Table 1 represent predictions of the peak pressure in 
reservoir 2* P 2 ma » on the basis of various viscometric assumptions,, In all 
cases the bulk modulus of the fluid* which is an important property in deter- 
mining P 2 © was obtained from measurements provided in Reference (2) for MIL-L- 
7808 fluid at 70°F o The bulk modulus was fit in the zero to 50*000 psi region 
with a third order polynomial given below* 

k = 2.739 x 10 5 + 1.827p + 7.996 x 10~ 5 p 2 - 7.615 x 10 _1 °p 3 

where p denotes the pressure in psi and k denotes the bulk modulus in psi* 

Column five of Table 1 shows values of P 2 max that would be predicted if the 
viscosity of the fluid were constant and equal to the measured value at 70°F 
of 23 centipoise. It is seen in all cases that when the viscosity is assumed 
to be constant at this value* the predicted values of P 2 max are consistently 
lower than the measured values* This result was somewhat surprising at the 
outset since it is well known that the viscosity of a. fluid increases with 
pressure and such increases in viscosity should result in more capillary re- 
sistance and still lower values of P2max° hence* the results shown in column 
5 indicate that if isothermal operation were occurring* the fluid would have 
to have a negative pressure coefficient of viscosity which is totally unrea- 
sonable* 

The above results suggested that thermal effects could indeed be influencing 
the results* The viscosity of a MIL-L-7808 oil was measured as a function of 
temperature over an interval of between 60°F and 300°F* The resultant (Fig* 31) 
viscosity-temperature curve was approximated with the viscosity index type 
formula given below* 

= exp { exp [ -3 * 7 05 £n (T + 460) 4- 24*39]} - 0*87* p - 14*7 

where the units for y e in the above expression are cp and T is in °F 0 

The above expression is accurate within the range of experimental accuracy of 1%, 
The next set of computations were carried out using the above viscosity temp- 
erature expressions 9 but still assuming no sensitivity to pressure* The results 
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are shown in column 6 of Table 1. Here in all cases the predicted pressure 
is substantially lower than the measured pressure, indicating that the true 
viscosity characteristics are somewhat higher than those assuming a zero 
pressure coefficient of viscosity. These results are somewhat more reasonable 
and considerably more encouraging. 

The final step in the analysis of the data shown in Table 1 was to assume an 
exponential pressure dependence on viscosity accompanied by the temperature 
dependence described above and to find the value of the pressure coefficient 
of viscosity, a, for which the predicted and measured values of P2 m ax were 
identical. Values of a so obtained which were assumed independent of temp- 
erature, are shown in the seventh column of Table 1. The values range from 
.47 x 10"^ to .75 x 10“^ in ^/lb with a mean value of .60 x 10~^ in ^/lb and a 
mean deviation of 0.114 x 10~^ in ^/lb. Values recorded in Reference (2) for 
a different batch of a similar fluid ranged from 10”^ in ^/lb at 68°F to 
.52 x 10 ^ in ^/lb at 210°F, approximately. Hence, values of the pressure 
coefficient of viscosity extracted from pulse measurements appear to lie with- 
in the range of those obtained from previously reported data. 

The above results indicate that the pulse viscometer appears to be a simple 
means for extracting pressure coefficients of viscosity. However, it is also 
seen that using reasonable values of pressure coefficients of viscosity the 
pulse data can be fit independently of any time delay phenomena and as such 
does not provide in itself a means for separating time effects from viscosity- 
pressure effects. In order to alleviate this difficulty, the dual capillary 
viscometer was constructed having two capillaries of different lengths. The 
different lengths result in significantly different residence times of a 
particle of fluid within the capillary and should provide relatively differ- 
ent delay behavior. This dual capillary viscometer arrangement with two 
secondary chambers (as shown schematically in Figure 32) which were connected by 
separate capillary tubes to the primary pressure chamber. The chamber volumes 
and capillaries were sized such that the geometry factor for both secondary 
chambers were nearly equal to one another. In this case, the diameter of each 
capillary was the same. The system provides two simultaneous pressure measure- 
ments. The first of these is the inlet pressure, pi(t), as a function of time 
which is the same for both capillaries. The second measurement is a pressure 
difference or the difference in values of p^(t) between the two exit reservoirs 
referred to as Ap 2 (t). 

In order to confine our attention to a situation providing as strong a time 
delay as possible, our efforts were concentrated on one pulse obtained with 
the use of a drop table at the pressure amplitude of the order of 20,000 psi. 

The pulse data as photographed on the oscilloscope is shown in Figure 33. In 
this case it is clearly seen that the input pulse is not sinusoidal, hence it 
was necessary to extract the data digitally from the photographs and feed in 
the input pulse, p^(t), into the computer to calculate the pressure difference 
Ap 2 (t). The two capillaries differ by approximately a factor of three in 
length and values of A p 2 ( t) correspond to values of P 2 (t) obtained from the 
short capillary subtracted from values of P 2 (t) obtained from the long capil- 
lary. The measured peak pressure difference, AP 2 max > is approximately 2800 psi. 

The first step in the analysis of this data was to compute the analytical 


C. 
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values of the pressure difference, Ap 2 (t) assuming no time delay, but includ- 
ing variations of bulk modulus with pressure and variations of viscosity with 
pressure and temperature . The results of these predicted measurements are 
shown in Figure 34. Measured and predicted values of ^P 2 max are plotted as 
functions of the pressure coefficient of viscosity. It can be seen that 
variations in the pressure coefficient of viscosity within a reasonable range 
cannot nearly account for differences between predicted and measured values 

of Ap 0 

2max 

At this stage, the full computer analysis described in the previous section 
was run with the inclusion of viscosity pressure, and temperature effects 
along with viscosity- time effects. The two triangular points shown in Figure 
34 correspond to calculated values of Ap assuming a constant delay time of 
t = 2 microseconds. 

time goes a long way toward reconciling predicted and measured values of 


Ap 


2max 


It is seen here that the inclusion of viscosity delay 
oward reconciling predicted and measured values of 
As calculations are costly here, no attempt was made to fully bracket 


the measurements of Ap 2rnax by varying the delay time. 


As a result of the temperature rises in the capillary^ it is possible that 
the viscosity can fall low enough to induce local turbulence. Such turbulence, 
if it did occur, could provide results in the same direction as those produced 
by short-time effects. It is doubtful, however, that turbulence could develop 
in the short times (less than 3 milliseconds) that the lubricant is subjected 
to the high pressure difference. The probability of turbulence would be 
greatly reduced with the use of smaller diameter capillaries and such experi- 
ments should be performed. The use of smaller diameter capillaries would also 
eliminate inertial effects which have not as yet been considered in the anal- 
ysis but could become significant at higher Reynold numbers. 
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SUMMARY AND CONCLUSIONS 


A pulse type capillary viscometer has been developed which appears to be 
capable of predicting pressure coefficients of viscosities with a simpler 
type of device than has heretofore been available. In addition, the apparatus” 
offers attractive prospects with regard to extracting viscosity-time character- 
istics. Preliminary experiments with a MIL-L-7808 oil indicate a delay time 
to exist in the order of 2 microseconds which is a qualitative range of the 
viscosity delay or time effects observed in ela stohydrodynamic rolling disc 
tractions . 

There are, however, several problems associated with the test apparatus. 

Thermal effects are prevalent as a result of viscous heating in the capillary, 
however they are believed to be fairly analyzable. There is also a possibil- 
ity of local turbulence occurring in some of the capillaries at some of the 
higher pressure differences. 

On the basis of the above conclusions, it is recommended that further modi- 
fications be performed to insure that there is no turbulence in the apparatus 
and to attempt to reduce thermal effects. One such method for doing this 
would be to apply an over pressure to the apparatus so that viscous character- 
istics at higher pressures could be investigated without resorting to extreme 
pressure differences. This would both reduce the thermal problem, inertial 
effects, and the tendency for turbulence. The use of an over pressure together 
with the ability to subject the apparatus to different ambient temperatures 
would provide a versatile viscometer for extracting more conventional viscosity- 
pressure characteristics under long times in a simpler manner than has been 
heretofore possible along with providing the ability to extract short time 
viscosity characteristics. 
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APPENDIX 1 


EQUIPMENT LIST 

One Harwood Engineering pressure intensifier Model E-1986 

Two Vishay Inter technology Incorporated strain gage bridges and 
amplifiers Model BAM-1 

One Tektronix Model 502A dual beam oscilloscope 
One Weston digital voltmeter Model 1420 
Micro -measurements foil gages #VM“SS-125CH-048 

Charles T. Gamble Incorporated subminiature resistors, manganin 
#200200AL, non-inductive wound 

One Spectra Physics helium-neon laser Model 132 

One Mechanical Technology Incorporated Fotonic sensor Model KD-38 
Machined Items per ten MTI drawings Series 264 
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NOMENCLATURE 


A 

a 

c 

E 

F 

f 

G 

G 


"1 

g(p) 

* 

g 

k 

L 

M 

P 

V P 2 


2 . 

2max 


t 

t 


u 

u 


v r v 2 


Cross-section of piston 

Input pulse amplitude 

Specific heat 

Mechanical pulse energy 

Viscosity- temperature function 

Dimensionless pressure in reservoir 1 

Dimensionless viscosity-pressure function 

Dimensionless equilibrium viscosity-pressure 
function 

Value of G in reservoir 1 

Dimensionless viscosity function 

Dimensionless viscosity function 

Bulk modulus 

Capillary length 

Mass of piston 

Pressure 

Pressure reservoirs 1 and 2 
Dimensionless pressure in reservoir 2 
Maximum value of p^ 

Flow rate 

Radius of capillary 
Radial coordinate 
Temperature 
Reference temperature 
Temperature in reservoir 1 
Time 

Delay time constant 
Pulse duration time 
Flow velocity 
Average flow velocity 
Volume reservoirs 1 and 2 
Axial coordinate 


. 2 
m 

lb/ in^ 

in- lb/ lb /°F 
m 

in- lb 
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NOMENCLATURE (Continued) 


a 

a 

P 

6 

6 


Ap 2 

Ap 2max 


M* 

M- 

M* 

M- 

P 

P 


e 

0 

1 


o 


T 


T 


d 


T 

max 



$ 

v 


0 


^ .c „ 

r I UbbUL C LOC1IXC1CUL U1 V-LO\_WO-A.i__y 

Dimensionless pressure coefficient of viscosity 
aa 

Dimensionless flow parameter 

Viscosity- temperature coefficient 

Dimensionless vise os it y- temperature 
coefficient, a6 

Pressure difference between discharge reservoir 
Maximum value of Ap^ 

Dimensionless capillary delay parameter 
Viscosity 

Equilibrium viscosity 

Reference viscosity at atmospheric pressure 

Viscosity in reservoir 1 

Density 

Reference density 

Dimensionless time D t/t Q 

Dimensionless delay constant tjt 

do 

Dimensionless time at which p^ is maximum 
Heat generation rate/unit length 
Heat generation rate/unit volume 
Dimensionless pressure function 


. 2 


/i k 


V 1 

psi 

psi 

lb- sec / in^ 

,lb-sec/ in^ 

lb- sec /in^ 

lb- sec /in^ 

lb /in 3 
m 3 

lb / in 
m 


lb- sec 
lb/ in^/sec 


28 



FIGURES 


29 




Tubing 




Schematic of Measured 


Primary Pressure Pulse 

Secondary Chamber Pressure 
Buildup 



mber Pressure Pulses 


31 


MTI-13628 




Secondary 

Chamber 

Housing 


Capillary 

Tubing 


Threaded 

Connector 


Coil Type Pressure 
Transducer Element 


Fig. 3 Secondary Chamber Components Showing Capillary Tubing 
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Gross Sectional View of Typical Secondary Chamber 
and Feed Through Capillary Tubing 













Digital Millivoltmeter 





Bridge and Oscilloscope Instrumentation 






Fig. 7 High Pressure Chambers 
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Fig. 8 Pulsed Pressure Data 

= 3,700 psi 

P 2 = 690 psi 

t = 2.65 milliseconds 
o 

Capillary: C 

Lubricant: 7808 

Run #28 



Fig. 9 Pulsed Pressure Data 

p^ = 7,000 psi 

p^ = 1,300 psi 

t = 2.35 milliseconds 
o 

Capillary: C 

Lubricant: 7808 

Run #29 
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Fig. 10 Pulsed Pressure Data 
= 10,400 psi 
P 2 = 1,970 psi 
t = 2.4 milliseconds 

Capillary: C 

Lubricant: 7808 

Run #32 


Fig. 11 Pulsed Pressure Data 
p^ = 14,750 psi 
P 2 = 4,950 psi 
t Q = 4.2 milliseconds 

Capillary: C 

Lubricant: 7808 

Run #40 
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Fig. 12 Pulsed Pressure Data 


Reproduced from 
best available copy. 


p^ = 7,600 psi 

p 2 = 840 psi 

t = 2.30 milliseconds 
o 

Capillary: C 

Lubricant: Jet II 

Run #58 


Fig. 13 Pulsed Pressure Data 

= 10,000 psi 

P 2 = 2,700 psi 

t = 5.50 milliseconds 
o 

Capillary: C 

Lubricant: Jet II 

Run #66 
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Fig. 14 Pulsed Pressure Data 


= 18,700 psi 

p^ = 4,000 psi 

t = 4.7 milliseconds 
o 

Capillary: C 

Lubricant: Jet II 

Run #67 


Fig. 15 Pulsed Pressure Data 
p^ = 31,000 psi 
p^ = 19,200 psi 
t Q = 6.6 milliseconds 

Capillary: E 

Lubricant: Jet II 

Run #84 
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Fig. 16 Pulsed Pressure Data 


= 8,400 psi 

P 2 = 1,170 psi 

t = 4.7 milliseconds 
o 

Capillary: E 

Lubricant: XRM 214B 

Run #95 



Fig. 17 Pulsed Pressure Data 

Pl = 9,200 psi 

P 2 = 1,450 psi 

t = 4.9 milliseconds 
o 

Capillary: E 

Lubricant: XRM 214B 

Run #90 
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Fig. 18 Pulsed Pressure Data 

= 22,200 psi 

P 2 = 8,500 psi 

t = 8.3 milliseconds 
o 

Capillary: E 

Lubricant: XRM 214B 

Run #100 


Fig. 19 Pulsed Pressure Data 
p^ = 27,000 psi 
P 2 = 10,300 psi 
t = 8.0 milliseconds 

Capillary: E 

Lubricant: XRM 214B 

Run #101 
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Fig. 20 Pulsed Pressure Data 


= 9,300 psi 

p^ = 6,000 psi 

t = 5.8 milliseconds 
o 

Capillary: B 

Lubricant : Polyphenyl 

Ether 

Run #135 



Fig. 21 Pulsed Pressure Data 

p^ = 12,200 psi 

p^ = 8,800 psi 

t = 5.4 milliseconds 
o 

Capillary: B 

Lubricant : Polyphenyl 

Ether 

Run #119 
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Fig. 23 The Variation of Peak Pressure 
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t -* 


Fig. 28 Schematic of Viscosity Response to a Step Jump 
in Pressure (t^ > % 2 > t^) 
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Viscosity, cp 




Impact Device 



• Left capillary longer than right capillary 

# V 1 >:> < V 2 + V ) 

Fig. 32 Differential Pulsed Viscometer Schematic 
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Fig. 33 Pulsed Pressure Data 

= 20,000 psi 

P 2 ~P 3 = 2,800 psi 

t = 3.8 milliseconds 
o 

Capillaries: D & E 

Run #113 
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